The Algebraic Bethe Ansatz and Tensor Networks 



O 



u 

■i— > 

a 

fl 
O 

o 



> 

in 
o 

<N 



V. Murg, 1 V. E. Korepin, 2 and F. Verstraete 1 

1 Vienna Center for Quantum Science and Technology, Faculty of Physics, University of Vienna, Vienna, Austria 
l C. N. Yang Institute for Theoretical Physics, State University of New York at Stony Brook, NY 11794-3840, USA 

(Dated: January 27, 2012) 

We describe the Algebraic Bethe Ansatz for the spin 1/2 XXX and XXZ Heisenberg chains with 
open and periodic boundary conditions in terms of tensor networks. These Bethe eigenstates have 
the structure of Matrix Product States with a conserved number of down-spins. The tensor network 
formulation suggestes possible extensions of the Algebraic Bethe Ansatz to two dimensions. 



I. INTRODUCTION 

The coordinate Bethe ansatz^ is an extremely suc- 
cessful method for solving one-dimensional problems ex- 
actly. It reduces the complex problem of diagonalizing 
the Hamiltonian to finding the solutions of a set of alge- 
braic equations. Once solutions to these algebraic equa- 
tions are found - numerical approaches to find them ef- 
ficiently exist in many cases - the eigenvalues are known 
exactly. However, the eigenstates are available only as a 
complex mathematical expressions the structure of which 
is not evident. This makes it insuperable, in general, to 
get interesting properties out of the states - like their en- 
tanglement characteristics or their correlations. The al- 
gebraic Bethe ansatz^ reveals more about the structure 
of the eigenstate and offers new perspectives to obtain 
scalar products^, norms^ and correlations^. 

In this paper, we point out this structure by formulat- 
ing the algebraic Bethe ansatz in the pictoresque tensor 
network language. In addition to making the ansatz more 
vivid, the tensor network formulation might bear the po- 
tential of extending the ansatz to higher dimensions. 

The description of states in terms of tensor networks 
has been very successful in the recent past. The one- 
dimensional matrix product states (MPS) 4,5 form the 
basis for the extremely successful density matrix nor- 
malization group (DMRG)£i Also, they have attracted 
considerable interest in the interdisciplinary field of 
quantum information and condensed matter physics^—. 
For describing the ground state of systems on higher- 
dimensional lattices, the projected entangled pair states 
(PEPS) 12 were introduced and proved to be useful for 
the numerical study of ground states of two-dimensional 
system o 13 i 14 . The Multiscale Entanglement Renormal- 
ization Ansatz (MERA) 15 i 16 allows the description and 
numerical study of critical systems. 

From the tensor network desription of the Bethe eigen- 
states it is immediately obvious that eigenstates can be 
described as MPS: see also Katsura and Maruyama [l7j . 
Katsura and Maruyama also show that the alterna- 
tive formulation of the Bethe Ansatz by Alcaraz and 
Lazo p~8l-[2Qj| is equivalent to the algebraic Bethe ansatz. 

In Sec. Ull we describe the tensor network form of the 
Bethe eigenstates and the structure of the obtained MPS. 
In Sec. [Till we formulate the algebraic Bethe ansatz in 
the tensor network language. In Sec lIV[ we give a pic- 
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FIG. 1. Tensor network constituting the Bethe eigenstate of 
the Heisenberg model or XXZ model with periodic boundary 
conditions. 



toresque description of the algebraic Bethe ansatz with 
open boundary conditions in terms of tensor networks. 



II. MATRIX PRODUCT STATE FORM OF 
BETHE SOLUTIONS 

Typically, Bethe-eigenstates are obtained as products 
of operators B(fij) applied on a certain vacuum state 
| vac ), i.e. 

|^(/ii,...,/i M )) = B(in)--'B(iiM)\vac). (1) 

The parameters {fij} are thereby solutions of Bethe equa- 
tions and the £?(/ij)'s play the role of creation opera- 
tors. In case of the antiferromagnetic Heisenberg model 
and the XXZ model with periodic boundary conditions, 
the vacuum corresponds to the state with all spins up 
and each operator B(jHj) creates one down-spin. Thus, 
the product of M such operators applied to the vacuum 
creates a state with M down-spins, i.e. magnetization 
S z = N/2 — M (with N being the number of spins). 
B(X) is an operator acting on the whole Hilbert-space of 
dimension 2^, but it has the well-structured form of a 
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Matrix Product Operator (MPO)2i with virtual dimen- 
sion 2. Indeed, as will be shown in Sec. IIII| 

s(A) = Yl (oi4 i (A)---4;wiiK 1 i ®---®of; 

Ii---In 

with fc,Z e {0,1}, of = |fc)(Z| (0 =t, 1 =i) and C k (X) 
being 2x2 matrices dependent on the parameter A. The 
product of operators B(/ii) • • • B(hm) can be read as the 
contraction of the set of 4-index tensors [C k {l^j)\ r r > with 
respect to a rectangular grid, as shown in Fig.[TJ Thereby, 
r, r', k and I label the left, right, up and down-indices, 
respectively. Explicitely, the matrices C k (A) read 

£ ° (A) = ( J C (°A) ) > £ ?( A ) = ( b (X) ) 
C . W = (°W) i£i(A)= (c(A)0) 

In case of the Heisenberg model Hxxx — EjLi ^xxx^ 
with 

^III = ^ [era; (8) cr^ + cr y (g) cr y + cr^ <g) o z - 1] , 
the functions 6(A) and c(A) are 

6(A) = — c(A) = — ^— . 
v ; 1 + A' v ; 1 + A 

For the XXZ model H X xz(A) = T,?=i ^ ] (A) with 

hxxz(A) = ^[a x ® a x + a y ® a y + A (a z ® cr z - I)} , 

the functions read 

= sinh(2zyy) = sinh(A) 

1 } sinh(A + 2zt7) ' 1 } sinh(A + 2ir]) ' 

The parameter rj is related to the inhomogenity A in the 
XXZ model via A = cos(2?7). 

Because of its "creation operator" -property, there is an 
inherent structure in the MPO B(/jl): each summand in 
the MPO B(/jl) must be non-zero only if k\ + • • • + fcjv = 
h + • • • + In + 1- This global constraint can be reduced 
to the local constraint that the tensors [C k must be 
non-zero only if r' = r + (fc — Z). This allows to inter- 
prete the virtual indices as "creation-annihilation" coun- 
ters: the right index r' is equal to the left index r if the 
physical state is unchanged, it is increased if a down-spin 
is created and is decreased if a down-spin is annihilated. 
Thus, the virtual indices transfer the information on how 
many down-spins are created and anniliated from left to 
right. Since the left boundary- state is (0 | and the right 
boundary-state is | 1 ) , it is guaranteed that the whole 
MPO creates exactly one down-spin. With the restric- 
tion r, r' G {0, 1} there are 6 possible configurations that 
fulfill the local constraint. In other words, only 6 entries 
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TABLE I. Disintegration of the matrices forming the MPS 
at step m, | ^ m ), into blocks. The full size of the matrices is 
D x D. The MPS has a conserved number of m down-spins. 



of the tensor [C k (/i)]£/ are non-zero. These 6 non-zero 
entries are 

[£8(m)]8, [c{(ri]\ 

[£S(m)]J, [^(m)]? 

[£g(M)]i, [£KM)]g, 

which is consistent with the matrices written above. 

The multiplication of all MPOs with the product 
state | vac ) evidently yields a Matrix Product State 
(MPS)&±£ with bond-dimension 2 M . Since each MPO 
B(X) has the "creation operator" -property to create one 
down-spin, the MPS contains exactly M down-spins. Ex- 
plicitly, the MPS reads 

|*) = J2 (0\(0\A kl ---A k "\0)\M)\k u ...,k N ) 

with matrices A k being block-diagonal in the sense that 
(a\(s\A k \P)\s') = [A k ]^ s sl . a and f3 are the virtual 
indices that range from to D — 1 (with D being the vir- 
tual dimension of the state). One the other hand, s and 
s' are the symmetry indices that transfer the information 
about the number of down-spins from left to right. The 
local constraint that guarantees this information trans- 
fer is s' = s + k. This constraint determines the blocks 
[wA^Ig, that are non-zero and allows a sparse storage 
of the state. The left boundary-state ( | and the right 
boundary-state | M ) fix the total number of down-spins 
of the MPS to M. 

The MPS is constructed iteratively by applying the 
MPOs B(ni), . . . , B(hm) successively to the vacuum 
state | vac ) . The state after m multiplications is evi- 
dently a MPS with m down-spins which shall be denoted 
as 

\*m)= £ (0|(0|^---^|0>|m>|fc 1 ,...,fc N > 

with A^ being block-diagonal in the sense that 
{a\{s\A k m \p)\s') = [A k m )f s , and fulfilling the constraint 
s' = s + fc, as before. The application of the operator 
B(/i) to | \I/ m ) yields a state with m + 1 down-spins 

l* m+1 >= (o |(o Kvi-^m+il 0)1^+1)1^,..,^). 

ki •••/cat 
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FIG. 2. Tensor network constituting the Bet he eigenstate 
of the Heisenberg model or XXZ model with open boundary 
conditions. 



The matrices *4^ + i emerge from tensor-products be- 



tween Cf and AL 



i.e. 



Ak 



= J2l£l®An- In indeX 

notation, 

(a\(r\(s\A k m+1 \f3)\r')\s') = Y / {r\C?\r')(a\(s\A l \f3] 

I 

Because of the constraints s' = s + 1 and r' = r + (k — £), 
5 = s + r and 5" = + r' suggest themselves as new 
symmetry indices. With this definition, S f = S + fc, as 
desired. 5 and 5" range from to m + 1, since 5 G 
{0, . . . , m} and r G {0, 1}. For S = and S = m + 1, 
there is the unique choice for 5 = r = and s = m, r = 1, 
respectively. For < 5 < m + 1, either 5 = 5, r = 0, or 
s = S — 1, r = 1. In this case, the index r must be kept to 
resolve this ambiguity. The index r can be incorporated 
into a new virtual index a as a = (a, r). Thus, the 
dimension of the blocks doubles for < S < m + 1. The 
column indices 5", s' ', r' and /? can be treated in the 



same way: for 5" = and 5' 



m 



1, s' and r' are 



unambiguously defined; for < S' < m + 1 there is an 
ambiguity that can to be resolved by incorporating index 
r' into a new virtual index j3 = (/?, r'). The matrices 
*A^ +1 in terms of the virtual indices a and f3 and the 
symmetry indices S and i.e. 

(a\(S\A k m+1 \P)\S'):=(a\(r\(s\A k m+1 \P)\r')\s'), 

have the desired block-form that fulfills the constraint 
S' = S + fc. Please refer to Table [J to see the dimensions 
of the block-representations that arise for different ra's. 

In the case of open boundary conditions, the Bethe 
Ansatz has the same form as in (pQ), merely the cre- 
ation Operators are not single MPOs, but products of 
two MPOs^^: 



B\- S (n) has the property to create 1 — s down-spins, 
whereas B s (ji) creates s down-spins (s G {0,1}), such 
that B(/i) is a creation operator for exactly one down- 
spin, as before. In terms of the previously defined 2x2 
matrices the MPOs read (see Sec. HV|) 

l\ ---In 

and 

s 1 _ s ( M )= ( s i4 1 M T ---4;(M) T io) f 1 i ®---® f; 

l± ■••In 

The virtual indices of B s {n) indicate the balance of cre- 
ated versus annihilated down-spins from left to right. 
This is due to the local constraint on [£f(/i)]JI/ that 
r' = r + (k — /), as mentioned before. Since the left 
boundary- vector is ( | and the right boundary- vector 
is | s ), the creation of s down-spins is guaranteed. In case 
of 5i_ s (/i), the MPO is built from the transposed matri- 
ces Cf (/i) T , such that the local constraint on [Cf (/i) T ]p is 
r = r' + (k— 1) and the virtual indices count the creation- 
annihilation balance from right to left. With the right 
boundary vector | s ) and the left boundary vector ( 1 | , 
^/one down-spin is created for s = and the number of 
down-spins is kept invariant for s = 1. 

The tensor-network representation for the Bethe-state 
with open boundary conditions is shown in Fig. [2j It con- 
tains twice as many rows as the tensor-network for peri- 
odic boundary conditions, which makes the contraction 
more challenging, in principle. However, as we see nu- 
merically, after a multiplication with a MPO-pair B(/i), 
the Schmidt-rank of the state only increases by a factor 
of 2 - not 4, as expected. This suggests that there should 
exist a representation with virtual dimension 2 also in 
the open boundary conditions-case. 



III. THE ALGEBRAIC BETHE ANSATZ 

Even though there exist numerous excellent reviews 
about the Algebraic Bethe Ansatz^^ - — 5 we resketch 
here the Ansatz in the picturesque Tensor- Network lan- 
guage for sake of completeness. In this way, it is trace- 
able, how the tensor networks shown in Figs. [1] and [2] 
form exact eigenstates of integrable systems. 



A. The Yang-Baxter Algebra 

In general, the starting point for the Algebraic Bethe 
Ansatz is the i?(A, /i)-tensor 



(2) 



s=0 



with a,f3,a',f3' ranging from 1 to some "auxiliary" di- 
mension d and A, \i being some complex parameters. This 
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FIG. 3. (a) Visualization of the 4-index i?-tensor R(\,/jl). 
In the abbreviated version, it is visualized as two crossing 
arrows with A attached to the up-down arrow and \i attached 
to the down-up arrow, (b) Yang-Baxter algebra as 4-index 
tensor with two virtual indices (left-right) and two physical 
indices (up-down), (c) Defining equation for the Yang-Baxter 
algebra. 
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FIG. 4. (a) Yang-Baxter equation in tensor- network form, 
(b) Abbreviated version. 



tensor defines the model under study, as will be shown 
later. Graphically, the tensor is represented by two cross- 
ing arrows, as shown in Fig. [3^, where A and ji are asso- 
ciated to the up-down and down- up arrows, respectively. 
After joining indices (af3) and (a' /?'), the tensor ([2]) can 
also be interpreted as matrix i?(A, /i) acting on the vector 
space V®V (with V = C d ). 

The condition on the i?-tensor (|2]) is that it fulfills 



FIG. 5. Inversion of the ordering of 3 composed Yang-Baxter 
algebras using i?-tensors. The inversion can be achieved in 
two ways, which makes necessary that the i?-tensors fulfill 
the Yang-Baxter equation (Fig. 



Yang-Baxter equation (star-triangle relation). Writing 

# (23) = 1 <g> R 
R^ = ii<g>l, 

the Yang-Baxter equation reads 

= R^ 2 \^ V )R^\\ V )R^ 12 \\ii). 

The graphical representation of this equation is shown 
in Fig. Another requirement is that solutions of the 
Yang-Baxter equation are regular, meaning that there 
exists a Aq and a z^o, such that 



#^3' (^o,^o ) - 8%'8p,- 



(3) 



The tensor i?(A, /i) defines the Yang-Baxter algebra 
T", (A) (a, a' = 1, . . . , d) by the relation 



R% (A, fiT*. (\)TL ( M ) = T£ ( M )T|, (X)R a J „ (A, M ) 

As usual, common indices are summed over. Defining the 
Monodromy T(A) as the matrix of operators 



Pi 



(Tl{\) ■■■ T](A)\ 



T(A) 



the definition of the Yang-Baxter algebra can be written 

as 

R(\, M ) [T(A)®T(m)] = [T( M )®T(A)] R(X, fi), (4) 

where the outer product "<g>" acts in the space V ®V 
in the sense that [T(/i)®T(A)]^, = T^{p)T^{\). T(A) 



FIG. 6. (a) Deformation of the Yang-Baxter equation, such 
that it yields the fundamental representation of the Yang- 
Baxter algebra, (b) Connection between the fundamental 
representation and the i?-tensor. (c) Formation of a more 
complex representation out of several fundamental represen- 
tations. 



can be considered as a 4-index tensor: 2 "virtual" indices 
a, a' of dimension d select the operator T^(fi) within the 
matrix, and two "physical" indices operate as input- and 
output index of the operator. T(A) is represented graph- 
ically in Fig. [3}d. The virtual indices are indicated as 
horizontal arrows; the physical input- and output indices 
are indicated as vertical in- and outgoing double- arrows. 
Using this graphical notation, the definition of the Yang- 
Baxter algebra assumes the simple form shown in Fig. [3k. 

In this picture, i?(A, /i) has the property to permute 
the thensors T(A) and T(/i). There is, however, one 
ambiguity that arises: there are two ways to go from 
T(A)<8>T(/x)<8>T(i/) to T(i/)<8>T(/z)<8>T(A). This inversion 
of the ordering can be achieved either by exchanging 
firstly A «-)• /i, secondly A o v and thirdly \i <H> z/, or 
by exchanging firstly v «-)• /i, secondly A O v and thirdly 
A o fi. This situation is depicted in Fig. [5] Thus, both 

(/i, u)R^ (A, u)R^ (A, /i) [T(A)®T(/z)®r(i/)] 
= [T(A)0T(/i)0T(z/)] ^ 12 ) (/i, l)^ 23 ) (A, ^ 12 ) (A, fj) 

and 

i?( 23 ) (A, /i)^ 12 ^ (A, ^)^ 23 ) (i/, /i) [T(A)(g)T(/i)0T(z/)] 
= [T(A)0T(/i)(g)T(z/)] i^ 23 ) (A, fx)fl< 12 ) (A, ^)i^ 23 ) (i/, /i) 

must be fulfilled. These two equations, however, are com- 
patible, because i?(A, fi) was required to fulfill the Yang- 
Baxter equation. This makes the definition of the alge- 
bra T^,(X) consistent. 

One representation of the Yang-Baxter algebra is easy 
to obtain - which is the fundamental representation. This 
representation is formed by the operators L^, (A, v) acting 



FIG. 7. (a) Co-multiplication property: formation of a new 
representation of the Yang- Baxter algebra out of two known 
representations, (b) Proof that the new representation still 
fulfills the defining equations for the Yang-Baxter algebra. 

on C d defined as 

[L%,(\,u)]f=R^ l (\,u). (5) 

In the graphical picture, the operators correspond to a 
clockwise "rotation" of the i^-tensor by 45 degrees, as 
shown in Fig. [6)3. The two indices attached to the hori- 
zontal arrow then become the virtual indices of the oper- 
ator, and the vertical arrow carries the physical indices. 
That these operators are a valid representation is due to 
the fact that the defining equation 

R(X, 11) [L(A, z/)(g)L(/i, v)\ = [L(fjL, z/)(8>L(A, v)\ R(\, /x) 

is just a "distortion" of the Yang-Baxter equation, as 
shown in Fig. [BJl. Up to now, the parameter v in L(A, v) 
is arbitrary. Most conveniently it is to set v = z/q. 

Once one representation L(X) is known, more com- 
plex representations are obtained by concatenating the 
L(A)'s horizontally, as depicted in Fig. [6k. Here, opera- 
tors T£>(A) acting on (C d )® N are constructed out of TV 
simple operators L^,(A) acting on C d via 

T£(A)= Yl L2 2 (\)®L%(\)®...®L a J(\). 

ck2,...,q;jv 

The outer product "(g)" affects the physical indices. In 
index notation, the operators read 

The operators defined in such a way fulfill (j4]), be- 
cause the i?-tensor subsequently interchanges the opera- 
tors L"/(A) from left to right - as can be retraced from 
Fig. [Tt for N = 2. Defining the matrices £f (A) as 
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FIG. 8. (a) Definition of the transfer matrix, (b) Proof that 
all transfer matrices commute. 

(a \Cf (A) | a' ) := [Z£,(A)]f , the operators T£,(A) assume 
the form of MPOs, 

W)= £ («l4 1 1 (A)---/:t Ar (A)k , K 1 ^---0of;, 



•/cat 



with of = \ k)(l\. 

The main building block of the Algebraic Bethe Ansatz 
is the transfer matrix t(A), obtained as the trace of the 
algebra T« (A), 

t(X) :=tr{T(A)}^^T-(A). 

The transfer matrix t(X) corresponds to T%,{X) with con- 
tracted left and right indices a and a' (see Fig. |5Jl). In 
the MPO picture, £(A) is represented by an MPO with pe- 
riodic boundary conditions. 28 Due to equation (j4]) that is 
fulfilled by the algebra, the transfer matrix has the prop- 
erty that [t(A),t(/x)] = for all A and /i. The way this 
property emerges from (J2J) can immediately be read off 
from Fig. [8]o: starting out with the expression £(A)£(/i), 
the identity in the form 1 = i?(A, /i) _1 i?(A, fi) can be in- 
serted at the virtual bonds; secondly, i?(A, fi) can be used 
to exchange T(A) and T(/i); thirdly, the cyclic property 
of the trace can be used to elimiate i?(A, fi) and i?(A, 
in order to end up with t(/x)t(A). 

This property makes £(A) the generator of an infinite 
set of commuting observables: if t(X) is Taylor-expanded 
with respect to A, t(X) = Io + Xh + A 2 i2 + • • then 
[/j, /jfe] = for all j and k. If one of the I^s is equal to 
the Hamiltonian of a model, it is called integrable, since 
there exist infinitely many symmetries which commute 
mutually. In fact, any function of t(X) can be generat- 
ing function for a set of commuting observables, like e.g. 
J-(X) = logt(A). The Taylor-expansion of this function 
reads 

T(X) = J-(Ao) + (A - Ao)J-'(Ao) + O ((A - A ) 2 ) . 



FIG. 9. (a) Precondition on the i?-tensor: at some point A = 
Ao, the i?-tensor decomposes into the outer product of two 
identities, (b) Logarithmic derivative of the transfer matrix 
t(A) at the point A = Ao- The first row represents t(Ao) -1 , 
the second row t'(Ao). 



It turns out that ^'(Aq) is local in the sense that 



^(A ) 



dX 



logt(A) 



N 



A=A 



+ 1) 



with only acting on sites i and i + 1. Thus, an 

integrable model is obtained described by a local Hamil- 
tonian 



TV 



ff = V>( 



M+l) 



Thereby, 



ft* 1 * 2 



A=A 



(6) 



(7) 



or 



h —R(\,v ) 



X=Xo 



respectively. To see this connection, it has to be realized 
that due to the regularity condition (j3j) t(Ao) is equal to 
the cyclic shift operator that shifts the whole lattice to 
the right by one site. The total momentum operator P 
is related to the cyclic shift operator according to 



:*(A ). 



(8) 



Graphically, t(Ao) is built from L", (Ao) shown in Fig. [9^i. 
The way the local Hamiltonian H emerges by differenti- 
ating the non-local expression .F(A) is sketched in Fig. [9b . 
Since •F'(Ao) = t(Ao) _1 t / (Ao), the first row in the figure 
corresponds to the inverted cyclic shift operator t(Ao) -1 



7 



and the second row corresponds to the derivative £'(Ao). 
The derivative £'(Ao) disintegrates into a sum of N deriva- 
tives with respect to each of the tensors L(A) at sites 
j = 1, . . . , N. As can be seen in the figure, term j has 
only support on two sites j and j+1 and thus corresponds 
to a two-site term that is related to the derivative L f (Xo) 
as formulated in ([7j). 

Models that emerge in such a way from combinations 
of fundamental representations are fundamental models. 
Examples are the spin- 1/2 Heisenberg model and XXZ 
model. In both cases, d = 2 and the .R-matrix assumes 
the form 



#(A,/i) 



(\ 



\ 



6(A,/i) c(A,/i) 
c(A,/i) 6(A,/i) 



\ 



1/ 



(9) 



Also, 6(A, /x) and c(A,/i) are of difference form, i.e. 
6(A, /i) = 6(A — /i) and c(A, /i) = c(A — /i). This yields a 
R- matrix of difference form, as well: i?(A, /i) = i?(A — /i). 
Explicitly, the functions b and c read 

b(X) 



c(A) 



1 + A 
A 



1 + A 



for the Heisenberg model and 



6(A) 
c(A) 



sinh(2i?7) 
sinh(A + 2i7?) 
sinh(A) 



(10) 

(11) 



sinh(A + 2ir]) 

for the XXZ model. Evidently, in both cases, R(0) = 
1, such that Ao = 0. In case of the Heisenberg model, 
#'(0) = hxxx with 

^iii = ^ [<r x (8) cr^ + cr y <g> cr y + cr^ ® a z - 1] . 

In the XXZ-case, 

hxxz(A) = ^[a x ® a x + a y ® a y + A (a z a z - t)} 

is obtained via R'(0) = 1/ s'mh(2ir])hxxz(A) with A = 
cos(2?7). 

Models (fundamental and non-fundamental) with R- 
matrix ([9|) are gl{2) generalized models. The Bethe 
ansatz for these models is especially simple and will be 
described in the following. 



B. Bethe Ansatz for gl(2) generalized models 

The Yang-Baxter Algebra with i?-matrix ([9]) is gen- 
erated by only 4 elements, such that the Monodromy 
assumes the form 



T(A) 



A(X) -B(A) 
C(X) D(X) 



with 

A(X) = T°(A), C(X) = Tq (A) 
B(X) = T?(A), -D(A) = Z?(A) ' 

The most important commutation relations of the alge- 
bra are 



B(X)B((i) = B(p)B(X) 
1 



A(\)B(jj) 



c(/x, A) 
1 

c(A,/i) 



B(M)A(A) - 4^B(A)A( M ) 



B(jx)D(X) 



c(m, A) 
c(A,//) 



The precondition for the Ansatz is that a representation 
must exist, for which there is a pseudo- vacuum |vac) 
that is an eigenstate of A(X) and D(X) and that is anni- 
hilated by C(A): 

A(A)| vac) = a(A)| vac) 
D(A)| vac) = d(A)| vac) 
C(A)|vac) =0. 

The goal is to diagonalize the transfer matrix t(X) = 
A(X) + D(X). Since all transfer matrices commute, 
[t(A),t(/i)] = 0, they have a common system of eigenvec- 
tors. Thus, all eigenvectors are independent of A. The 
eigenvalue problem reads 

t(A)|*)=r(A)|*>. 

The Bethe Ansatz 

|*(^i,.,.,/xm)> = B(m)---B(n M )\vac). 

fulfills the eigenvalue problem provided that the /i/~'s ful- 
fill the Bethe equations 



d(fJLn) 



M 

n 



(12) 



(n=l,. . . ,M). The eigenvalue r(A) is then equal to 

M M 

The proof is obtained by utilizing algebraic relations only 
and can be gathered from appendix lAl From r(A), the 
eigenvalue of the Hamiltonian (j6j) is obtained as 



E 



t'(Ao) 
r(A ) 



(13) 



The total momentum is, according to (jSJ), equal to 

p=-ilnr(A ). (14) 
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C. Bethe Ansatz for the Heisenberg model and the 
XXZ model 



The addends are usually referred to as magnon momenta 
that can be written as 



In case of the Heisenberg model and XXZ model, this 
spezializes as follows: the matrices Cf (A) that build up 
the MPOs T^/(A) have block form. Written out, they 
read 



£g(A) = 



££(A) 



1 

c(A) 

b(X) 





, £?(A) 




c(A) 
1 



These MPOs are symmetry conserving in the sense that 
T£, (A) changes the number of down-spins by a' — a. This 
is due to the local constraint that [Cf (A)]^, are non-zero 
only if a' = a + (k — Z), as discussed in Sec. [TTJ 

Using these considerations, the vacuum state is obvi- 
ously the state with no down-spins, namely 



| vac) = |0) 



^|0) 



(0 =t, 1 =|) • This state is annihilated by C(A) and is 
an eigenvector of A (A) and D(X) with eigenvalues 



o(A) = 1, 
The Bethe- Ansatz state 

| *(//!,. ,.,/xm)) 



TV 



d(A) = c(A) 



B(m)--'B(n M )\vac) 



is a state with M down-spins, i.e. with magnetization in 
z-direction equal to S z = \N — M. 

The Bethe equations obtained by the Algebraic Bethe 
Ansatz are equal to the equations obtained by the Coor- 
dinate Bethe Ansatz. In case of Heisenberg model, it is 
advantageous to introduce variables Zj that are related 
to /ij in (p2j) via 

Zj 1 

^ = 2i " 2 

for a direct comparison with results of coordinate Bethe 
ansatzi22r— : In terms of these variables, the Bethe equa- 
tions read 



N 



M 

n 



2% 



2i 



(15) 



with n = 1, . . . , M . From the Bethe solutions {zj}, the 
energy is obtained using ([T3|) as 



E = 



r ; (0) 
r(0) 



M , 



According to (fH|) . the total momentum P has eigenvalue 

M 



P 



-ilnr(0) = ( -^ ln 



2 7 " + z 



p j = 7r — 2 arctan(zj). 



using the identity 

arctan(,) = ilni±^. 

In term of the magnon momenta pj , the total momentum 
reads 



M 



(16) 



and the energy is equal to 

M 

£=-2^(l-co S fc)). 
i=i 

For solving the Bethe equations (fT5j) it is advantageous 
to bring them to their their logarithmic form 

M 

Np n = 2nl n + ^ e(p n ,pj), 

i=i 

where 

, 0(p,g) » a 

2 cot = cot cot - . 

2 2 2 

and are integers G {0, . . . , N}. Solutions can then be 
found iteratively, as described in (3l|. The ground state 
configuration for N even and M = N/2 is (ii, . . . , Jm) = 
(1,3,...,A/--1). 

In case of XXZ model, it is advantageous to introduce 
the variables Zj related to jij in ([T2|) via 

H = z 3 -iV + i^ 

to compare with the Coordinate Bethe Ansat z 33 i 34 . The 
Bethe equations then read 



cosh(z n — ir])\ N yr sinh(z n — Zj — 2irj) 
cosh(z n + if] J -fjj sinh(z n — Zj + 2ir\) 



From the Bethe solutions {zj}, the energy is obtained as 



sin(2?7) 2 



(0) cos(2t?) + cosh(2^) " 

The total momentum obtained from ([T4|) is again of 
form ([16]) with 

^ = —2 arctan (tanh(zj) tan(^)) . 
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(a) 

2 A 

(b) x-? 
(d) 
(e) 



A 3 





P(A) 



P(A) 



FIG. 10. (a) Permutation symmetry, (b) Yang-Baxter equa- 
tion, (b) Partial transpostion symmetry, (e) Crossing unitary 
condition, (d) Unitary condition. 



In terms of the momenta pj , the energy can be expressed 
as 

M 

£ = -2^(A-cos( Pj )). 
The Bethe equations in their logarithmic form read 

M 
3 = 1 



that permutes two indices. Using the matrix- notation 
i?(A, /i) from appendix IIII1 i.e. considering the .R-tensor 
as matrix acting on V V with V = C d , a variant of 
the i?-tensor with the first two indices permuted can be 
defined as 

n(x,fjL) = pr(\h) 

The basic assumption is that the i?-tensor fullfills the 
symmetry condition 

PR(\n)P = R(\li) 

(see Fig. ITOk). Then, the i?-tensor can expressed just by 
two crossing arrows and it is not necessary to distinguish 
between them by marking them with the arguments. In 
fact, it is assumed in the following that R is of difference 
form, i.e. i?(A, /i) = R(X — /i). Thus, the tensor R(X — 
fi) will be characterized by two crossing arrows together 
with the argument A — /i, as shown by the rightmost 
depiction in Fig. [TOk . Using this notation, the Yang- 
Baxter equation assumes the form shown in Fig. [TUb . 

It is furthermore useful to define the partial transposi- 
tion 



which is equivalent to flipping the direction of "up-down" 
arrow. In analogy, 



corresponds to flipping the direction of the "down-up" 
arrow. Accordingly, the partial transposition symmetry 
condition 



with 



cot 



0(p,q) 



Asin^ 



cos2±2 - Acos^ 



(17) 



and Ij G {0, . . . , N}. The ground state configuration for 
N even and M = N/2 is again found with (Ji , . . . , Im) = 
(1,3,...,7V-1). 



IV. ALGEBRAIC BETHE ANSATZ FOR OPEN 
BOUNDARY CONDITIONS 

The method described for periodic boundary condi- 
tions is generalizeable to models with open boundary con- 
ditions and boundary fiek k) 22 i 23 i 35 . We resketch here the 
Ansatz for open boundary conditions following closely 
Sklyanin 23 using a picturesque language. 

For the following it is required that the i?-tensor fulfills 
several conditions. To express these, it is convenient to 
define the permutation operator 



P = Y,\3,i){i,3 



ft(A)* 1 = Tl(X) t2 

is expressed by Fig. ITUb . 

Further conditions are the unitarity condition 

n(X)TZ(-X) = p(X) 

and the crossing unitarity condition 

^(A) tl ^(-A-2c) tl = p{X) 



(18) 



(19) 



1,3 



with p(A) and p(A) being some scalar functions of A and 
c denoting some constant characterizing the .R-tensor. 
These conditions are represented by Figs. fTUll andflQfe. 



A. Reflection Algebras 

As the Bethe Ansatz for periodic boundary conditions 
is based on the Yang-Baxter algebra, the footing of the 
open boundary conditions Ansatz are the reflection al- 
gebras JC~(X) and /C + (A) spanned by {/C~£(A)|a, \3 = 
1, . . . , d} and {/C + (A) Q! ^|a, j3 = 1, . . . , d}. The graphical 
representation of these two algebras is shown in Fig. [TTh : 
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(a) 



(b) 



4 



A 



. A - /,, 



—A - n - 2c j 

—X + u y y 



1L 



(a) r(A)=tr(X:-(A)^+(A)) 

(b) 



r(A)r( M ) 



It 



"A^u /-A + 



' 1J 



A 


A - M / -A + yj 

* /A + M / 

/ -A - At - 2c 















4 ifc-'ii 



: T (/i)r(A) 



FIG. 11. (a) Graphical representation of the reflection alge- 
bras AC - (A) and /C + (A). The horizontal arrows indicate the 
virtual indices a and f3; the vertical arrow indicates physical 
indices, i.e. the input- and output indices of the operators 
/C^(A) and /C^g(A) respectively, (b) Defining equations for 
the reflection algebras (reflection equations). Each intersec- 
tion of two lines represents an i?-tensor. The argument of the 
i?-tensor is written next to the intersection. 



as in the case of the Yang-Baxter algebra, each of the 
two algebras is considered as a 4-index tensor with 2 
"virtual" indices a and j3 of dimension <i, represented 
by the horizontal arrows, and the 2 "physical" indices 
(corresponding to the input- and output indices of the 
operators JC~a(X) and [/C^g(A)] respectively), represented 
by the vertical arrows. The only difference to the Yang- 
Baxter case is that the virtual indices both are on the 
right-hand side of the tensor in case of K~ (A) and on the 
left-hand side in case of /C + (A). The correspondence to 
the Monodromy in the open boundary condition case is 
the matrix of operators 



/C±(A) 



//C±(A) ••• /C±(A)\ 



U±(A) ••• /C±(A)/ 



The denning equations for the reflection algebras are 
the reflection equations, represented by the tensor net- 
work in Fig. [TTb . In this figure, each intersection of two 
lines represents an i^-tensor. The argument of the R- 
tensor is written next to the intersection. Algebraically, 
the reflection equations read 



FIG. 12. (a) Definition of the transfer matrix, (b) Proof of 
the commuting property of the transfer matrix, [r(A), r(fi)] = 
0. 



and 

1Z(-X + /i)[/C+(A)] £l ft(-A -fi- 2c)[£+(/i)] t2 



with 

^(A) =/C ± (A)0l 

A^(A) = l(g)/C ± (A). 

The outer product "(g)" is thereby interpreted as in ((4]) 
and t is the d x d identity matrix. 

Using these algebras, it is possible to define a commut- 
ing set of transfer matrices via 

r(A) = tr (/C"(A)/C + (A)) . 

Graphically, r(A) corresponds to JC~ (X) and /C + (A) be- 
ing glued together, as shown in Fig. [T2r . The commu- 
tativity of the transfer matrices, [r(A),r(/i)] = 0, can 
be proven using the unitary and crossing unitary con- 
ditions ([18]) and (fT9|) and the reflection equations. The 
proof is sketched in Fig.[T2b: starting out with t(A)t(/x), 
the line connecting /C + (/i) and JC~ (fi) can be pulled over 
the line lying above that connects /C~(A) and /C + (A) us- 
ing (fT9]) and over the topmost arrow connecting the two 
A- algebras using JT5J). Next, the network is mirrored 
vertically by using the reflection equations. Finally, the 
drawn out line is pushed back using (fT8|) and (fT9|) , which 
leads to r(/x)r(A), as desired. 

What is remaining is to find concrete representa- 
tions of the reflection algebras. Examples of simple 
representations with physical dimension 1 have already 
been found^ More complex representations can be con- 
structed by assembling a known representation with two 
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(a) 



(b) 



B 

*■ 



(d) t(a) 



XL 

A 



(c) 



* 



A AAA A A A A 



A AAA A A A A 



▲ A 

A A 



4 A A 



AAA 

A A 



> A 



■- 

I I" 



A-,. 



A - 



.It 



7 a+m (a) a V 



/A — n 



(b) 



1> 






i r 






A 


A + u 






A + p 


(d) 

A — 


M 


> 




A 




A 


u 




A' 


u 


a-m\ 





FIG. 13. (a) Composition of a new representation of AC - (A) 
out of one pair of i?-tensors and a known representation of 
AC - (A) (that already fulfills the reflection equations). The 
known representation is indicated by the shaded surface, (b) 
Simple representation of /C + (A) with physical dimension 1. 
(c) Composition of a complex representations of AC - (A) of 
dimension d N by attaching N pairs of i?-tensors to a simple 
representation with physical dimension 1. (d) Transfer matrix 
built from the representations (b) and (c). 



FIG. 14. Proof that the composed representation shown 
in Fig. 113a fulfills the reflection equations: the main idea 
of the first two steps (a) and (b) is to pull the vertical line 
rightmost by applying the Yang-Baxter equation twice. The 
three i2-tensors to which the Yang-Baxter equation is applied 
are marked by the shaded triangles. The new situation now 
allows the application of the reflection equations, as shown in 
step (c). The last step (d) consists in pushing the vertical line 
back by applying the Yang-Baxter equation twice, such as in 
steps (a) and (b), but in reverse order. 



i^-tensors in the way shown in Fig. [T3k . The physical 
dimension of the new representation is thereby increased 
by a factor d. That this assembly is indeed a valid repre- 
sentation can be proven using the Yang-Baxter equation 
and the reflection equations. The proof is sketched in 

Fig.CE! 

Thus, starting out with a simple representation with 
physical dimension 1 for /C~(A), a representation with 
physical dimension d N is obtained after N iterations with 
the relation expressed in Fig. IT5k . The structure of the 
representation after N iterations can be gathered from 
Fig.[T3b. Assuming a simple representation with physical 
dimension 1 for /C + (A) (depicted in Fig.[T3b). the transfer 
matrix assumes the form shown in Fig. [T3l i. 

For the sake of simplicity, we choose the simple repre- 
sentations with physical dimension 1 equal to the identity 
(which is a valid representation that fulfills the reflec- 
tion equations). Using the notation for the fundamental 
representation of the Yang-Baxter algebra introduced in 
equation (|5j) and Fig. [6b, the representations of the alge- 
bras Kr (A) and /C + (A) look as shown in Fig. [T5k . The 
transfer matrix assumes the form depicted in Fig. [T5b . 
Algebraically, the representation of /C~(A) is then the 
product of two MPOs, 



K P W = Y. R MK^). (20) 

s=l 

In terms of the previously defined matrices Cf (/i), the 



(a) 



L ± 



r t 



(b) 



T T 
t t 




-L JL X 



r(A) 



T J 



r x i 



T T 



T T T 

N -2 N-l N 



FIG. 15. (a) Representation of /C~(A) and /C + (A) for open 
boundary conditions, (b) Transfer matrix built from these 
representations. 



MPOs read 



K 



(.i< i (M)---< Ar (M)i/^)< 1 



ki • • -/ctv 
h •••In 
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T 1 
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1 t t f 



FIG. 16. Derivation of the open boundary condition Hamil- 
tonian by derivative of the transfer matrix shown in Fig. 115b 
at the point Ao- (part I). 



t f I 



An 



! I 



T U 

1 



I X 









a; 

t 








t 



FIG. 17. Derivation of the open boundary condition Hamil- 
tonian by derivative of the transfer matrix shown in Fig. 115b 
at the point Aq. (part II). 



and 

km= E (^i4 1 (^) T ---4;(M) T i«)o>---®of; 

k±---kN 
Ii---In 

with of = | k )(l |. The representation of /C + (A) has phys- 
ical dimension one and is equal to the identity with re- 
spect to the virtual indices, i.e. 

^apW = 

The transfer matrix constructed in this way is indeed 
related to a local Hamiltonian with open boundary con- 
ditions. This Hamiltonian is obtained as the derivative of 



the transfer matrix at the point Ao at which the i?-tensor 
is equal to the identity (see (j3j)). Explicitly, the obtained 
Hamiltonian is of the form 



N-l 
i~l 

and related to the transfer matrix via 

1 



(21) 



H 



2d 



r'(Ao). 



In ([2T]h the symbol a refers to an auxiliary system that 
is traced out. Using the notation from appendix IIII] this 
relation is seen as follows: the derivative t'(Ao) disin- 
tegrates into a sum of 2N terms, each term containing 
one tensor differentiated at Ao and 2N — 1 tensors eval- 
uated at Ao- Due to the regularity condition (j3]) of the 
i?-tensor, the tensors evaluated at Ao assume the simple 
form shown in Fig. [9^. As can be gathered from Figs. [16] 
and[T7J each differentiated tensor at site i corresponds to 
a two-site term for i = 1, . . . , TV — 1 (with h being 

defined in Q). For i = TV, two indices of the tensor are 
traced out, which leads to the one-site term tr a h^ N,a \ 



B. Bethe Ansatz for the XXZ model with open 
boundary conditions 

For the XXZ model, the virtual dimension d is equal 
to 2, such that the Monodromy can be written in the 
form 



jc-(\) 



A{\) B{\) 
C(A) V(\) 



with 



.4(A) = /C - (A), C(A) 
6(A) = /C 1 (A), P(A) 



Ku(A). 



The i?-tensor has the form (j9j) with b(X) and c(A) being 
defined by ([TO]) and JTTJ). It fulfills the regularity condi- 
tion ([3]) at the point Ao = 0, the unitarity condition ([T8]) 
with p(A) = 1 and the crossing unitarity condition ([T9]) 
with c = 2ir] and p(X) = 1 — sin(2^) 2 / sin(2^ — i\) 2 . Us- 
ing representation ([20]) for /C~(A), the i?-tensor generates 
the Hamiltonian 



H 



1 



(H° x bc xz (A)-A) 



sinh(2i^) 



with 



H° x bc xz (A) 



N-l 

E< 

n=l 



IXXZ 



(A). 



The precondition for the Bethe Ansatz is that a rep- 
resentation must exist, for which there is a pseudo- 
vacuum | vac ) that is an eigenstate of A(X) and V{X) 
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and that is annihiliated by C(A): 

w4(A)| vac) = a(X)\ vac) 
V(X)\vac) = d(X)\vac) 
C(X)\vac) = 0. 

As argumented before, the operator C(A) annihilates one 
down-spin, whereas *4(A) and V(X) keep the number of 
down-spins constant, such that the state with all spins up 
is a valid pseudo- vacuum. The goal is now to diagonalize 
the transfer matrix r(A) = A(X)+V(X). Since all transfer 
matrices commute, [r(A),r(/i)] = 0, all eigenvectors are 
independent of A. The eigenvalue problem reads 

r(A)|*)=r(A)|*>. 

The Bethe Ansatz 

| ^(/ii, . . . ,/i M ) ) = • • -B{hm)\ vac). 

fulfills the eigenvalue problem provided that the /i/s ful- 
fill the Bethe equations. The proof is based upon the 
algebraic relations between *A(A), 6(A), C(X) and V(X) 
and is described in detail in [231 ] . 
Defining the momenta pj via 



the Bethe equations in their logarithmic form rea d 36 i 37 
{N+l)p n = 7r/ n +B(p n , -p n )+ — 

with Q(p,q) being defined in (fT7j) . The ground state for 
TV even and M = N/2 corresponds to (7i,...,/m) = 
(1,3, ... ,7V- 1). The energy eigenvalue of Hg% z (A) for 
a configuration (pi, . . . ,Pm) is obtained as 

M 

^ z (A) = -2^(A-cosfe)). 

V. CONCLUSIONS 

Summing up, we have sketched the Algebraic Bethe 
Ansatz using the pictoresque language of tensor net- 
works. In a future paper, the method will be extended to 
[three-dimensional] space lattices and its physical impli- 
cations for cohesion, ferromagnetism and electrical con- 
ductivity will be derived^ 
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Appendix A: Algebraic Derivation of the Bethe 
Equations 

For completeness, we sketch here the derivation of the 
Bethe Equations using algebraic relations. We thereby 
follow Korepin 2 . 

The goal is to find eigenvectors of t(X) = A(X) + D(X) 
using algebraic relations between A(X), 5(A), C(A) and 
D(X). The commutation relations that are required are 

B(X)B(ji) =B(ji)B(X) (Al) 
A{\)B{») =f(\,riB(»)A(\)+g(\,n)B(\)A(rtA2) 
D(\)B(n) =f( ft ,\)B( Jt )D{\) +g(v,\)B(\)D(^A3) 

Here, the abbreviations /(A,/i) = l/c(/i, A) and g(A,/i) = 
— 6(/i, A)/c(/i, A) are used. 
The Bethe Ansatz reads 

|*(^i,...,/x M )> = B(/n)---B(/i M )\vac), (A4) 

where | vac ) is a state that is an eigenvector of A(X) and 
D(X) with eigenvalues a(A) and d(A), and that is anni- 
hilated by C(X). A(X) applied to | ^(/ii, . . . , /lm) ) using 
relation (|A2|) yields in principle 2 M terms, because each 
commutation of A(X) with a B(/ik) yields 2 terms and it 
takes M commutations to move A (A) from left to right. 
However, these two terms are not arbitrary. Both terms 
only perform exchange operators: the /-term in (|A2[) ex- 
changes the operators A and £>, but not their arguments; 
the g-term, on the other hand, exchanges the operators 
A and B and their arguments. Due to this, after M 
commutations the following conditions must hold: 

• Every term must contain M £>'s and one A. 

• The M + 1 coefficients (A, /xi, . . . , /am) are dis- 
tributed among the M £>'s and the one A. 

Since all £?'s commute, there are only 2 cases: either A is 
argument of A - then the term looks like 

B(/i 1 )---B(/i M )A(X)\vac). (A5) 

Or, A is argument of one of the £>'s. Then the term is of 
the form 

B(X)Y[B(/i J )A(/i n )\vac) (A6) 

with n e {1,...,M}. Thus, the 2 M terms can be col- 
lected into M + 1 linearly independent terms: 

A(X)\ ^(/ii, . . . , ijl m ) ) = A B(/n) • • • B(/i M )A(X)\vac) 

M 

+ 5> n £(A) J] B(/i 3 )A(/i n )\vac) 

71=1 j'/:n 

What remains to be done is the calculation of the coeffi- 
cients A and A n . 
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The expression (|A5|) is obviously obtained after M 
commutations using the /-term in (|A2|) . The g-term 
must not be applied, because it introduces a B(X). Thus 

M 

To obtain (|A6|) . it is convenient to rewrite the Bethe 
Ansatz (fA4|) as 

|^(/ii,...,/i M )) =B(n n ) Y[B(nj)\vac). 

This is possible for all n, since all B's commute. Since 
expression (jA6|) must not contain £?(/i n ), the first com- 
mutation with A(A) must be performed using the g-term 
in ()A2|) . The expression then reads 

g(\/i n )B{X)A{/i n ) \\B(/ij)\vac). 

All further commutations must use the /-term, because 
another use of the g-term would introduce B(/i n ) in the 
expression again. Thus, the coefficients must be 

A n = g(X, /l n ) Yl /(/in, Mj) 

The application of -D(A) to | ^(/ii, . . . , /lm) ) can be 
treated in a similar way using relations (|A1|) and ((A3]). 
Again, the application yields M + 1 terms 

D(A)| . . . , /i M ) ) = A • • • 5(/i M )^(A)|Mc> 

M 

+ ^A n £(A) B^)D(/i n )|t;oc) 

n=l j=^n. 



The coefficients are 

M 

and 

An = #0n, A) ]^[ /(/ij,/i n ). 

Thus, I ^(/ii, . . . , //m) ) is an eigenvector of t(X) = 
A(X) + D(A) if 

a(/i n )A n + d(/i n )A n = 

for n = 1, . . . , M. These relations are the Bethe Ansatz 
equations, which can be written in the form 

d(Lln) TT c(/i n? /ij) 
a 0n) c(/ij,/i n ) 

under the assumption that g(X, /a) is an odd function in 
the sense that g(A, /i) = —g(/i, X) (as it is the case for the 
Heisenberg model and the XXZ model). 
The eigenvalue r(A) is obtained as 

r(A) =a(A)A + d(A)A, 

which can be expressed as 

M 1 M 
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